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We demonstrate, using first-principles calculations, that the electronic structure of FeSei-ccTe^ 
(a;=0.5) is topologically non-trivial, characterized by an odd Z 2 invariant and Dirac cone type 
surface states, in sharp contrast to the end member FeSe (x=0). This topological state is induced 
by the enhanced three-dimensionality and spin-orbit coupling due to Te substitution (compared to 
FeSe), characterized by a band inversion at the Z point of the Brillouin zone, which is confirmed 
by our ARPES measurements. The results suggest that the surface of FeSeo.sTeo.s may support a 
non-trivial superconducting channel in proximity to the bulk. 

PACS numbers: 73.43.-f, 74.70.Xa, 73.20.-r 


Among the Fe-based superconductors, the FeSei-^Te^ 
family of compounds m is of particular interest. First, 
it has the simplest PbO structure (space group PA/nmm) 
with Se (or Te) atoms forming distorted tetrahedra 
around Fe (see Fig.JlJa)) similar to the structure of FeAs 
planes in the families of FeAs-based high T c supercon¬ 
ductors [5 . Second, the internal parameters can be sys¬ 
tematically tuned by the substitution of Se by Te [SHE], 
which provides us a platform for in-depth study of pos¬ 
sible superconducting mechanisms and topological char¬ 
acters. Thirdly, superconductivity has been observed for 
a wide range of composition x mu, and the transition 
temperature T c can be further enhanced by pressure It 
[TTj . More recently, superconductivity with T c higher 
than 77 K was suggested for single unit cell FeSe films m 
epitaxially grown on SrTiOs substrates. 

Despite these interesting properties though, the par¬ 
ticularities of the system have still not been fully ex¬ 
plored. Earlier studies, both theoretical and experimen¬ 
tal, suggest the similarity of the electronic structures of 
the Fe chalcogenides (FeSe, FeTe) p~3lfl5] and the FeAs- 
based [46lfT8] superconductors. Indeed, the low-energy 
physics around the Fermi level is dominated by the Fe-3d 
states, and the morphology of the Fermi surfaces is sim¬ 
ilar. On the other hand, a surprisingly stable (no split¬ 
ting under external magnetic field) zero-energy bound 
state (ZBS) at randomly distributed interstitial excess 
Fe sites was observed in very recent scanning tunneling 
microscopy (STM) measurements on the surface of su¬ 
perconducting Fe(Te,Se) [19], suggesting possible topo¬ 
logical feature of its electronic structure. Obviously, the 
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5 p orbitals of Te are more extended and have stronger 
spin-orbit coupling (SOC) than the 4 p orbitals of Se. The 
consequences of Te substitution, particularly for the bulk 
topological character of FeSei-^Te^, have been largely 
ignored in the literature and will be the main purpose 
of the present paper. Based on first-principles calcula¬ 
tions combined with angle resolved photoemission spec¬ 
troscopy (ARPES) measurements, here we report that 
the electronic structure of FeSeo.sTeo.s is topologically 
non-trivial, in sharp contrast to its end member FeSe. 
The topological properties of FeSeo.sTeo.s can be charac¬ 
terized by an odd Z 2 number, and the existence of Dirac 
cone type surface states, in proximity to bulk super¬ 
conductivity, should favor topologically superconducting 
surface states, as suggested by Fu and Kane m ■ 

The electronic structures of FeSe and FeSeo.sTeo.s 
are calculated with SOC included. The calculations 
are performed based on the density functional theory 
(DFT) [2T1I22] and the generalized gradient approxima¬ 
tion (GGA) for the exchange-correlation potential [23] . 
as implemented in the plane-wave pseudopotential based 
BSTATE (Beijing Simulation Tool of Atomic TEchnol- 
ogy) package [24]. The experimental lattice parame¬ 
ters [6] are used in the calculations. Maximally localized 
Wannier functions (MLWF) [25] are constructed from 
bulk calculations, and then used to study the surface 
states in the semi-infinite system. To treat the substi¬ 
tution properly, we have calculated FeSeo.sTeo.s by us¬ 
ing both the virtual crystal approximation and the two- 
formula cell with ordered Se and Te sites. Both calcula¬ 
tions give converging results. 

We first neglect the SOC and concentrate on the com¬ 
parison between the electronic structures of FeSe and 
FeSeo.sTeo.s (as shown in Fig. [l] (c) and (d)). The band 
structure of FeSe is very similar to that of LaOFeAs as 
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FIG. 1: (Color online) The crystal structure and DFT elec¬ 
tronic band structures of FeX (X = Seo.sTeo.s)- (a) is the 
crystal structure with the :c-axis pointing along Fe nearest 
neighbors. Grey and green balls represent Fe and X atoms, 
respectively, (b) Schematic plot of the hybridization along the 
z axis between the combined orbitals D~ y and P~ (see details 
in appendix B), consisting of intralayer pda and interlayer ppn 
hopping, (c) Band structure of FeSe with internal parameter 
zx=0.2345 without SOC. (d) Band structure of FeSeo.sTeo.s 
(zx=0.2719) without SOC. (e) Band structure of FeSeo.sTeo.s 
with SOC. (f) Zoom-in view of the solid red box area in (e). 
The size of red the circles in (c) and (d) indicates the weight 
of the p z component of the chalcogen atoms. The two d z 2 
bands are indicated below Ep in (d). The original A 5 bands 
split into A.q and A 7 . Along the Y -Z line, two A 6 bands cross 
and hybridize to open a SOC gap of about 10 meV. The red 
dashed line corresponds to a Fermi curve across the gap. 


reported before mi- At the T point, the valence band top 
is not occupied, leading to the well known hole pockets of 
Fermi surfaces around Y. The three top-most states at T 
can be labeled as two-fold degenerate Y$ states ( d yz /d xz 
orbitals) and non degenerate Y^ state ( d xy orbital), re¬ 
spectively. There is a clear band gap larger than 0.6 
eV above the valence band top at the Y point. All the 
d — d anti-bonding states (with negative parity) are lo¬ 
cated above the gap except the d z 2 orbital that has a 
weaker anti-bonding state. Among them, the most in¬ 
teresting state is the second highest one with remarkable 
red circles, which belongs to the representation and 
comes from the anti-bonding d xy orbitals of Fe and p z 
orbital of chalcogen. The weight of the p z component is 


illustrated by the size of red circles, which suggests that 
the state can be affected by hybridizing with the d xy 
and p z orbitals. Looking along the Y — Z direction, the 
band dispersion of the T ^ state is the strongest among 
all d states. Nevertheless, since it is energetically high, 
the band structure around the Z point of the Brillouin 
zone (BZ) is only slightly affected and remains similar 
to that around Y. The result suggests that FeSe is quite 
two-dimensional with strong anisotropy. 

Immediate differences can be seen in comparing these 
results with the band structure of FeSeo sTeo 5 shown in 
Fig# 1 ): (1) the r 2 state in FeSeo.sTeo.5 is significantly 
pushed down and almost touches the valence band top, 
and the band gap at Y (above the valence and top) is 
nearly closed; (2) the band dispersion of this state 
along the Y — Z direction is strongly enhanced. As a 
result, a band inversion occurs at the Z point, which 
implies a change of topological property. 

Compared to FeSe, the c-axis and the experimental 
intra-layer distance d z of FeSeo.5Teo.5 are enlarged by 
8 .2% and 11 %, respectively (See appendix A), while the 
a-axis and the inter-layer distance change only very lit¬ 
tle. As a result, the intra-layer hybridizations, espe¬ 
cially the pdcr hopping as shown in Fig. [ljb), are seri¬ 
ously weakened, and the band center is lowered very 
much and becomes close to as shown in Fig. [ljd). 
On the other hand, because the Te 5 p orbitals are much 
more extended than the Se 4 p orbitals, the inter-layer hy¬ 
bridization through the ppir bonding becomes stronger. 
Therefore, the Te substitution enhances the hopping be¬ 
tween layers and gives rise to a larger dispersion for in 
FeSeo.5Teo.5- These two observations are in good agree¬ 
ment with the band structure of FeSeo.5Teo.5 shown in 
Fig. 0 (d). Therefore, the band inversion in FeSeo.sTeo.5 
is a consequence of the weakened intra-layer hopping 
and enhanced inter-layer hopping originating from the 
Te substitution. 

Once SOC is included, the band structure of 
FeSeo.sTeo.5 opens a direct SOC gap and a non-trivial 
Z2 invariant can be defined by assuming a “curved chem¬ 
ical potential”- the red dashed line in Fig. 0 e)-lying 
between the 10th and 11th bands (neglecting the spin- 
doublet degeneracy of the bands). Generally, the dou¬ 
bly degenerate Y 5" states split into Yq and Y ^, and the 
odd state turns into Yq (see details in appendix B). 
Along the Y-Z high-symmetry line, two A@ bands un¬ 
der symmetry hybridize and open a gap of about 10 
meV, which is clearly shown in Fig. 0t)- When defin¬ 
ing a Fermi curve through the SOC gap, the Z2 invari¬ 
ant is easily calculated from the parity criterion, which 
comes out to 1 . (The parities at all time-reversal invari¬ 
ant momenta (TRIM) are presented in appendix B.) This 
non-zero Z2 invariant indicates that FeSeo.sTeo.s is in a 
topological phase that can support the non-trivial surface 
states (SS). Due to the substantial SOC of Te, increasing 
the content x enlarges the SOC gap, which is benefi- 
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FIG. 2: (Color online) ARPES spectra along the T-M direc¬ 
tion in the a polarization geometry, (a) FeSeo.sTeo.s before 
surface evaporation, (b) FeSeo.sTeo.s after 6 minutes of K 
surface doping, (c), (d) MDCs corresponding to data from 
the dash boxes in panels (a) and (b), respectively. 


FIG. 3: (Color online) Normal emission ARPES spectra. 
First row: schematic band dispersions. Second row: ARPES 
data on FeSeo. 5 Teo. 5 - Third row: ARPES data on LiFeAs. 
First column: dispersion at k x = 0 along the T-Z direc¬ 
tion. Second, third and fourth columns: electronic dispersions 
along cuts indicated in the first column. 


cial for the detection of the SS in the gap. Moreover, 
we also performed dynamical mean-field theory (DMFT) 
calculations to confirm the band inversion and identify 
the strong band renormalization due to electronic cor¬ 
relations (see appendix C for details). The correlation 
effects do not change the detail of the electronic bands, 
but simply reduce the bandwidth. 

We performed ARPES measurements in order to 
demonstrate experimentally the existence of the band 
with strong p z orbital character crossing Ep along the T- 
Z direction. Large single-crystals of FeSeo.45Teo.55 were 
grown using the self-flux method, and LiFeAs with FeAs 
flux method. ARPES measurements were performed at 
the Advanced Light Source and Synchrotron Radiation 
Center, using a VG-Scienta electron analyzer. The light 
used was linearly polarized in directions parallel to the 
analyzer slit. The K source used for evaporation is made 
of a SAES K dispenser. In the experiments, the largest 
coverage is less than one mono-layer. All data were 
recorded with linear horizontal polarized photons with 
a vertical analyzer slit (a geometry). Under this config¬ 
uration, odd orbitals with respect to the emission plane 
are visible, such as the d xy and d yz orbitals, while the d xz 
band should not be detected. In addition, our experimen¬ 
tal setup leads to a z component of the light polarization, 
and thus orbitals with mainly z components such as d z 2 
and p z orbitals are also observable [ 26 ]. 

The comparison of the DFT band calculations on FeSe 
and FeSeo.sTeo.s, shown in Figs. [ljc) and 0d), suggests 
that while it is pushed far above Ep at T in FeSe, the 
band in FeSeo.5Teo.5 forms an electron band just above 


Ep. To prove its existence, we raised the chemical po¬ 
tential by performing in-situ K doping, and we measured 
ARPES spectra along T-M at hu = 30 eV, which coin¬ 
cides with k z = 0 m- The results before and after evap¬ 
oration (6 minutes) are shown in Figs.[ 2 ja) and[ 2 jb), re¬ 
spectively. As expected, the hole band is sinking further 
below Ep after evaporation. Interestingly, an additional 
electron band is observed, as clearly shown by contrast¬ 
ing the momentum distribution curves (MDCs) obtained 
before (Fig.[ 2 jc)) and after (Fig.[ 2 jd)) evaporation. This 
band, locating about 30 meV above the top of the va¬ 
lence band in FeSeo.sTeo.s, is very similar to the small 
3 D electron pocket reported in (Tl,Rb) 2/ Fe2- a; Se2 [ 26 ] . 
which mainly has a p z component. 

In the presence of SOC, the DFT calculations indicate 
that the Tg”, Tj and Tj states at the T point mostly 
come from the d xz , d yz and d xy orbitals, respectively, 
while the Tq state (labeled as without SOC) have 
an important p z component besides the d xy orbital. As 
shown in Fig. [ 3 ja), while the p z band locates above Ep 
at T, it is shifted below Ep upon moving along T -Z and 
it reaches its minimum at the Z point. On its way down, 
the p z band crosses the d xy , d yz and d xz bands, opening 
a SOC gap with the d xz band. This situation is illus¬ 
trated schematically in Fig. [ 3 ](a) . We now ask how the 
hybridization of the p z band with the d xz band should 
affect the ARPES measurements. Around the V point, 
the p z band forms an electron-like band located above 
Ep along the T-M direction. Therefore, it should not 
be observed in our experimental geometry. Due to selec- 
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tion rules, the d xz band should not be observed either 
and only the d xy and d yz band can possibly be seen, as 
shown schematically in Fig. | 3 |b). However, because it 
hybridizes with the p z band, one should expect to be 
able to detect the d xz band near the crossing point be¬ 
low Ep, as illustrated schematically in Fig.[ 3 jc). Finally, 
our calculations indicate that the p z band has a hole-like 
dispersion near the Z point, and we should be able to 
observe it. On the other hand, the d xz band no longer 
hybridizes strongly with the p z band and its intensity 
should be significantly suppressed again in the a polar¬ 
ization, as described in Fig. [ 3 ^d) . 

In order to confirm this dispersion, we compare in 
Fig.^e) ARPES spectra recorded on FeSeo.45Teo.55 with 
different photon energies. We see some intensity between 
the T and Z points that we assign to the p z band sinking 
down from T to Z. At the Z point, this band has merged 
with the strong d z 2 band and it is thus undistinguishable. 
We display in Figs. [ 3 jf)-[ 3 jh) three ARPES intensity cuts 
along T-M recorded on FeTeo.55Seo.45 and corresponding 
respectively to k z values around T (Cut#l, Fig. [ 3 ^f )), 
the hybridization between the p z and d xz bands (Cut#2, 
Fig. j 3 ^g)), and around Z (Cut# 3 , Fig. [ 3 jh)). As the 
d xy band heavily renormalizes compared to the d xz /d yz 
bands, it is shallower and hardly resolved with weaker 
intensity [ 28 ] . In Cut#l (Fig. [ 3 jf )), as expected, none 
of the p z and d xz bands are detected near k z = 0 . In 
contrast, not only the d yz band is observed in Cut#2 
(Fig.^g)), but the d xz band as well due to its hybridiza¬ 
tion with the p z band. At the Z point, away from hy¬ 
bridization, the d xz band disappears from the ARPES 
spectrum (Fig. [ 3 ](h)). Unfortunately, the p z band is too 
close to the strong d z 2 band to be distinguished unam¬ 
biguously. 

As shown in Figs.[ 3 ji) to [ 3 ^ 1 ), the strong k z dispersion 
of the p z band is also observed in LiFeAs, which exhibits 
very clear spectral features. Because its bottom is located 
away from the d z 2 band, the dispersion of the p z band 
along T-Z can be identified very clearly, as illustrated in 
Fig. | 3 |i). As with FeSeo.45Teo.55, the intensity of the d xz 
band is the strongest along Cut ^5 (Fig. | 3 [k)), where it 
hybridizes with the p z band. Figs. [ 3 ]^k)^[ 3 jl) also show 
clearly that the p z band has a hole-like dispersion below 
Ep’t as expected theoretically. 

Finally, we analyze the spin-resolved Fermi surfaces 
around the Dirac point of the semi-infinite system formed 
by the SS to identify the non-trivial topology. With the 
surface Green’s function calculated from the modified ef¬ 
fective Hamiltonian considering DMFT on-site modifica¬ 
tion, the spin-filter surface states and the corresponding 
Fermi surfaces can be obtained directly. From the dis¬ 
persion of the SS shown in Fig. | 4 ja), the protected SS 
emerges in the SOC gap around Ep due to the non-trivial 
topological nature of the bulk system. The Fermi surface 
of the SS (Ep = 50 meV) is illustrated as a bright yel¬ 
low circle in Fig. [ 4 ^b), and the spin orientation for the 




k x (2jt/a) 

FIG. 4: (Color online) Topological non-trivial 001 surface 
states and spin-resolved Fermi surface of these states, (a) Sur¬ 
face LDOS on the (001) surface for FeSeo.sTeo.s considering 
DMFT on-site modification, (b) Fermi surface of the topo¬ 
logical state (bright yellow circle beside the projected bulk 
states). The in-plane spin orientation is indicated by green 
arrows. 


SS around the Fermi surface is marked by green arrows, 
which form a n Berry phase enclosed. The magnitude of 
the z spin component is very small compared to the in¬ 
plane component. The 7r Berry phase signifies the topo¬ 
logical non-trivial properties of the bulk. By inducing a 
s -wave superconducting gap, the chiral SS can play an 
important role to produce the Majorana zero mode. 

In conclusion, we have presented both theoretical and 
experimental evidences for a topological non-trivial phase 
in FeSeo.sTeo.s. From the DFT calculations, we show 
that the band topology is sensitive to the intra-layer and 
inter-layer hopping terms, which can be tuned by the Te 
substitution. The Te substitution content x strongly af¬ 
fects the structural, electronic and topological properties 
of these materials. We have identified the topologically 
non-trivial electronic band structure of FeSeo.sTeo.s with 
a band inversion, characterized by an odd Z2 invariant 
and spin-moment locked SS. Our ARPES data strongly 
support that the band forms a band inversion at the 
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Z point. The similar results can also be applied to iron- 
pnictides such as LiFeAs. Due to the topological non¬ 
trivial surface states, the FeSei-^Te^ materials would be 
an ideal system for realizing possible topological super¬ 
conductor and Majorana fermions on the surface. 
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Appendix A: Structural parameters of FeX and band structure of LiFeAs 


TABLE I: Structural parameters of PbO-structure FeX. The lattice parameters are from experimental data [6], and both 
optimized (Opt.) and experimental (Exp.) internal chalcogen positions zx are shown. d z is the cartesian distance in the z 
direction from the X plane to the Fe plane. 


a (A) c (A) 

zx(Opt.)/d^ (A) 

zx(Exp.)/d z (A) 

FeSe 3.784 5.502 

FeTe 3.8215 6.2695 

FeSeo.sTeo.s 3.7933 5.9552 

0.2345/1.2902 

0.2496/1.5648 

0.2476/1.4745 

0.2652/1.4591 

0.2829/1.7736 

0.2719/1.6192 


For LiFeAs, we used the experimental lattice parameters [29] and relaxed the free internal coordinates. The 
experimental structure (space group P4/nmm, 129) has the Li sites in 2c positions, which lie above and below the 
centers of the Fe squares opposite the As. The calculated DFT band structure is given in Fig. [HJ The T^ band shows 
a large dispersion and a band reversion along the z direction. 



FIG. 5: (Color online) Electronic band structure of LiFeAs. The p z component is highlighted by the size of the red circles. 


Appendix B: Symmetry and parity analysis 

Since the system has inversion symmetry, it is convenient to combine these orbitals into bonding and antibonding 
states with definite parity as 

I Dt) = d=(\Fe a ) ±\Fe' a )), 

I Pp) = A(|^}t|^)) 

Using the orbitals with definite parity defined in the main text, we can label all the non-spin-orbital (NSO) bands 
with the irreducible representations (IRs), which are given in Table [nj The T^ band is composed of the antibonding 
states with both d xy and p z characters. At both the T and Z points, the IRs of the group are labeled as T^, 
whereas along the T-Z high-symmetry line, we label these IRs as A n because of its C± v symmetry. In our conventions, 
the (x,y) axes are rotated by 45° as compared to the crystallographic axes, so that the d xy orbital in our definition 
is the one pointing from Fe to chalcogen atoms, as shown in Fig. 1. 
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TABLE II: Combined orbitals and related irreducible representations 


D/P (+/-) 

z 2 

yz/xz 

xy 

2 2 

x -y 

x/y 

z 

D4 h + 

rf 

r t 

rt 

rt 

rt 

r t 

- 

r 3 

r 5 

r 2 “ 

rr 

rr 

r 2 “ 

+ 

A x 

A5 

A4 

A3 

A5 

Ai 

- 

a 4 

A5 

Ax 

A 2 

A5 

Ai 



FIG. 6: (Color online) Electronic band structure along T-Z. (a) The IRs are labeled in the NSO bands, (b) The parities are 
labeled in the SOC bands. 


The IR labels along the T-Z line are shown in Fig. j^a). Around Ep, the D+ z ^ yz bands can mix with the P+/ y 

orbitals with the same T^ representations, although the P-D hybridization is not strong due to the little orbital 
overlap in the real space. For the r^“ representation, the D~ character can mix with the P~ band, and the d xy and 
p z orbitals have a strong hybridization along the z direction. However, d z becomes larger with the Te substitution 
content x increasing, thus leading to weaker hybridization. As a result, the band, consisting of d xy and p z orbitals, 
sinks down to Ep and has a large dispersion in FeSeo. 5 Teo. 5 - Moreover, due to selection rules, the bonding state of 

Z line. As a consequence, the Tj band characterized by the 

Beside, all 

the IRs in NSO bands and parities in SOC bands are presented in Table |III| denotes that the quadruple 

bands consist of two Kramers pairs with opposite parity. The non-trivial Z 2 index is presented beside the vertical 
line, which denotes the hypothetical Fermi level in the main text. 


D+ y cannot mix with the p z band, including along the T 

D+ y representation exhibits a weak dispersion in the 2 direction, which is not sensitive to the height d 


Appendix C: LDA+DMFT confirmation 

The LDA + DMFT method has proven to be a powerful technique to study the electronic structure of correlated 
systems. In this section, we apply this method to FeSeo.sTeo.s with a local Coulomb integral U=Fo=4.0 eV and 
a Hund’s coupling J=0.7 eV, and confirm that it has a topological character of band inversion by computing the 
topological invariants within the DMFT framework. 

The one-electon spectral function is defined as 

1 ImE(u;) 

kCJ 7T (uj + /1 — ek — ReYi(uj)) 2 + /mS(w) 2 
























TABLE III: The IRs and parities for each band. 


IRs(NSO) 


p 



D 



r 

r t r 2 

r 5 r 5 


r t r t r s 

P 

r+ 

1 4 

F 17 r 5 - 

z 

r 2 - r+ r 5 - r+ r+ 

rf r 3 - 

r 2 “ 

r+ r+ 

1 4 1 5 

rr r 5 - 

Parities(SO) 

1 2 

3 4 5 

6 

7 8 9 

10 

11 12 

13 14 15 16 

r 

+ - 

- - + 

+ 

+ + - 

+ 

+ + 

- - 

z 

- + 

- - + 

+ 

+ + - 

- 

+ + 

+ - - - 

X 

(+ -) 

(+ -) (+ 

-) 

(+ -) (+ 

-) 

(+ -) 

(+ -) (+ -) 

M 

(+ -) 

(+ -) (+ 

-) 

(+ -) (+ 

-) 

(+ -) 

(+ -) (+ -) 

R 

(+ -) 

(+ -) (+ 

-) 

(+ -) (+ 

-) 

(+ -) 

(+ -) (+ -) 

A 

(+ -) 

(+ -) (+ 

-) 

(+ -) (+ 

-) 

(+ -) 

(+ -) (+ -) 

z 2 





1 



in terms of the LDA band dispersion ek and the self-energy E(cj). This momentum-resolved spectra Ak(cj) is shown in 
the upper panel of Fig. [7| where the overall renormalization of the bands and the bandwidth reduction are apparent. 
However, the relative positions of different bands almost do not change. The T^~ and T^ bands crossing near Ep 
remains in the correlated system. 


0.4 - 


■ 


- 0.4 - 




FIG. 7: (Color online) (Upper panel) LDA+DMFT calculations. (Lower panel) SOC band structure only taking the 
LDA+DMFT on-site modification into consideration. 

In a previous work |30|, only the Green’s function at zero frequency G^ x (0) = fi — e k — £k(0) is needed to 
determine the topology of the quasiparticle states, since renormalization does not change the non-trivial topological 
nature of the bulk system. Following Refs. EH. we compute the topological invariant of an interacting system with 
the Hamiltonian defined by H t (k,u;) = iL(k, cj) + £k(0) — /i. A modified effective dispersion of H t ( k) for which we 
neglected the imaginary part of £k(0) is given in the lower panel of Fig. [ 7 ] By introducing correlations, the bandwidth 
is reduced by approximately half, but the band inversion character still exists. According to the parity criteria, we 
conclude that FeSeo.sTeo.s remains a topological correlated system. The topological surface state is derived from this 
modified effective TB Hamiltonian. 















